Shapes and textures of ferromagnetic liquid droplets 



O 

o 



Shubho Banerjee* and M. Widom 
Department of Physics, Carnegie Mellon University, Ptttsburgh, Pennsylvania 15213 

(February 1, 2008) 

Theoretical calculations, computer simulations and experiments indicate the possible existence 
of a ferromagnetic liquid state. Should such a state exist, demagnetization effects would force a 
nontrivial magnetization texture governed by the shape of the liquid droplet. Since liquid droplets 
are deformable, the droplet shape couples to the magnetization texture. This paper solves the 
joint shape/texture problem subject to the assumption of cylindrical droplet symmetry. The shape 
undergoes a change in topology from spherical to toroidal as exchange energy grows or surface 
tension decreases. 
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■ I. INTRODUCTION 

(N 

In a spontaneously magnetized liquid state, long range magnetic order would exist in the liquid without application 
tj^ of any external field. The existence of a ferromagnetic liquid state in dipolar fluids has been indicated by mean field 
O \ calculations |^-|^ and computer simulations |l7|-|l0|l. Experiments to observe liquid ferromagnetism in ferrofluids pd| ] 
^ . are challenging because the fluids often freeze |12yi3| ] or phase separate well above the predicted low temperatures 



o . 

a ■ 



' for the onset of spontaneous magnetization. Recent experiments on strongly interacting FeaN ferrofluids |15| do show 
. a hint of a possible ferromagnetic transition. Similarly, experiments on super-cooled Co-Pd alloys indicate the 
possibility of ferromagnetism in supercooled liquid metals. In this case it is the strong exchange interaction, and 
'"^ not the dipole interaction, that would cause the spontaneous magnetization. The experimental evidence in both 

■ ferrofluid and supercooled Co-Pd is inconclusive, since the bulk of evidence addresses the temperature dependence of 
O \ paramagnetic susceptibility above the Curie temperature ]l5|-[l^ . Definitive proof of ferromagnetism below the Curie 

I ^ I > temperature remains elusive. 

Although the existence of a ferromagnetic liquid state is yet to be proven experimentally, spontaneous polarization 
. coupled with other order parameters has already been observed. Some electrically polarized liquid crystals |l9| ] 
^ ' show a helical ordering of the dipole moments in the liquid. In superfluid ^He the magnetic moment couples to the 
, superconducting order parameter |^^. Many superfluid "^He phases are therefore also magnetically ordered. 

■ It is interesting to consider the magnetization texture (spatial variation of the orientation of magnetization) inside a 
droplet of such a ferromagnetic liquid ||2l| . The magnetization texture likes to avoid poles ||2^ to minimize its energy. 
However, this leads to defects inside the texture. For example, the rotating magnetization texture with cylindrical 
symmetry inside a sphere 



M(r,0,0) = M^, (1) 

y where (j) is the unit vector for the variable (j), avoids all poles but has a vortex line running through the center. Near the 
' , \ vortex of such a texture the magnetization is topologically unstable and might escape into the third dimension p^ ] 
^ with a nonzero component along the vortex line. Whether this happens depends on the balance between demagnetizing 
Q and vortex energies. Simulated annealing of the magnetization inside a cubic box suggests that replacing vortices with 
^ , point defects may be favorable However, for sufficiently large droplets the demagnetizing energy will dominate 
^ ' and only textures with vanishing demagnetizing energy (per unit volume) will occur. 
• ^ Any defect is likely to have a system-shape-dependent energy cost causing a deformable liquid droplet to deviate 

from a spherical shape. The complete calculation of the shape of an unconfined ferromagnetic liquid droplet in three 
M J dimensions, coupled with the calculation of its magnetization texture, remains an interesting and challenging unsolved 
■ - - 1 problem. 

This problem has a simple solution in two dimensions in zero field. The magnetization texture inside any soft (zero 
anisotropy) ferromagnetic solid thin film is given by van den Berg's algorithm p5| which avoids all poles, and thus 
all magnetostatic energy, at the expense of a domain wall through the film. A liquid droplet, which can change its 
shape, prefers a circular shape to minimize its surface energy. The magnetization lines inside a circle form concentric 
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circles according to van den Berg's algorithm. For a circular shape the domain wall energy is also minimized because 
the domain wall shrinks to a point vortex. The circle thus solves the coupled texture and shape problem in zero field. 

In a previous paper |^ addressing the thin film limit, we described the evolution of magnetization texture and 
droplet shape under the application of an in-plane magnetic field. The vortex stretches to become a domain wall, which 
displaces towards the edge of the droplet. If exchange energy is taken into account, the droplet exhibits reflection 
symmetry-breaking, immediately revealing its magnetized state. 

Our present goal is to analyze the shape and texture of an unconfined ferromagnetic liquid droplet. The problem 
is highly nontrivial even in the absence of applied field. To simplify the study, we restrict our attention to shapes 
and textures with an axis of continuous rotational symmetry. We find the texture exhibits a planar character with a 
vortex line along the symmetry axis. When surface tension is high or magnetization weak, the droplet shape remains 
nearly spherical, with slight "bulging at the waist" and dimpling where the vortex line meets the surface (Fig. |^). 
As surface tension drops or magnetization grows, the dimpling increases, shortening the vortex line. Eventually the 
droplet undergoes a change in topology to a toroidal shape (Fig. ^). 



II. SHAPE, TEXTURE AND ENERGY 

We consider droplets of volume V and shape 17 with an internal magnetization texture M(r). Four terms contribute 
to the energy. 

The first term is the energy of the droplet surface. Oil, 

Esurf (fr^aA (3) 

Jon 

with A the total droplet surface area and a the surface tension which we take to be isotropic. The surface energy is 
minimized for a spherical droplet. 

The second term is the exchange energy, which, for an isotropic medium, can be written as the integral over the 
droplet volume of 

1 dMkdMk 

Uexch = T^O^^ ^ , (4) 

2 axi oxi 

where a is the exchange constant and the summation convention is employed [^7| . The exchange constant has units 
of length squared. For a metallic ferromagnet the exchange constant typically takes a value of order Xa^, where a 
is a characteristic atomic size and A ~ 10^ is the ratio of exchange field to magnetostatic field For a ferrofluid, 
we can derive an effective exchange constant reflecting the energy cost of placing finite size particles into a rotating 
magnetic texture. By calculating the correction to the mean field in a cavity of diameter a caused by rotation of 
the magnetization, we find a = Ana^/W. Again, a is of order Aa^, where now A = 1 because the energy cost is 
magnetostatic in origin. 

The third term is the energy associated with unmagnetized regions within the fluid. We presume that throughout 
most of the droplet the magnetization M(r) has constant magnitude, AI, that minimizes the free energy density f{M). 
Inspecting eq. ^ we see that exchange energy density may diverge at the core of a vortex, where Uexch ^ a/r'^. 
Within a distance Vc of the center of the vortex the exchange energy density Uexch matches the free energy density 
cost Ucore ^ /(O) — f(.M)- The core radius so-obtained does not depend on droplet shape or magnetization texture. 

The final term in (0) is the demagnetizing energy that arises as a consequence of the long range 1 /r^ character of 
the dipole interaction. This weak fall-off of the interaction makes the total energy of the system dependent on global 
magnetization texture and system shape in general, challenging our notions of thermodynamic limits ||2^ . The shape 
dependent demagnetizing energy for dipolar systems can be best understood in terms of the demagnetizing field 

Hz,(r) = / d'r' (M(r') • fi(r'))^-^ - / d'r' (V • M(r'))^-^. (5) 
Js |r-r'|'^ Jy Ir-r'l'^ 

Here fi is the normal to the surface of the system. The first term on the right hand side reflects the surface poles that 
are created where the magnetization has a component normal to the surface. The second term is the contribution 
from the bulk charge density that is created by a nonzero divergence of the magnetization. 

Hd is called the demagnetizing field because it inhibits the magnetization, as shown by the magnetostatic energy 
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ED = -\j cPvliD{v)-M{v), (6) 



which can be rewritten 



Ed = ^ fdh |Hz,(r) 

J all space 



(7) 



This energy is manifestly positive definite and hence inhibits the magnetization. 

To avoid a demagnetizing energy by ehminating its demagnetizing field, magnetization textures in crystals break 
up into domains separated by domain walls. For suitable domain structure the demagnetizing energy can be removed 
entirely. The domain wall width is set by balancing the cost in exchange energy for rotating textures against the 
cost in magneto-crystalline anisotropy energy when a the magnetization does not point along an easy axis. Because 
a ferromagnetic liquid should lack magneto-crystalline anisotropy, the domain wall width diverges pi] ] unless it is 
limited by some other characteristic length ||2^ . 

The four terms in eq. (|^) grow at differing rates as droplet volume increases while droplet shape and texture are 
held fixed. We find 

E,urf ~ crRh (8) 
Eexch ^ aM'^hlogR/rc 

Ecore ^ Ucoreh 
Edemag ~ DM^E?h. 

In the above, we have taken magnetization texture (|^) in a cylinder of height h and radius R for our calculation of 
Eexchi and we have taken M(r) constant inside an ellipsoid of demagnetization factor D for our calculation oi Edemag- 
For sufficiently large droplets, Edemag dominates unless a texture is found to reduce it or remove it altogether. 
Anticipating that textures with vanishing Edemag will emerge, the dominant energy becomes Egurf, so large droplets 
will favor compact shapes. For large R, we find Eexch dominates Ecore- At fixed i?, the exchange energy mimics the 
core energy, with an energy cost proportional to the vortex length. 

III. ENERGY MINIMIZATION 

We now confront the problem of simultaneously solving for the lowest energy magnetization texture within a droplet 
of a given shape while varying the shape to achieve the lowest total energy. To simplify our work, we choose to work 
within a limited subset of the family of possible textures and shapes. We impose cylindrical symmetry, motivated both 
by our suspicion that the true energy minimum may exhibit such a symmetry and by the considerable calculational 
simplifications that it allows. 

Hence we assume the shape Q. and the magnetization texture M(r) are symmetric under rotations about the z axis. 
The magnetization vector field obeys 




M(r,0,z)= sin0 cos0 M(r,0,z) (9) 



The droplet shape is defined by its boundary dVL which we parameterize by the function ±z(r). This family of 
shapes exhibits reflection symmetry through the z = plane in addition to rotational symmetry about the z axis. 

For large droplets the demagnetizing energy is dominant, as we showed in eq. (^. Consider the problem of 
minimizing Ejj within a symmetric shape. Since Eo > 0, we can solve this minimization with any texture that 
achieves Ed = 0, which, by equation requires Hp — through all space. Inspecting we can achieve H^) = 
by eliminating all surface poles and volume divergence. 

The simplest texture with H/j = is given in eq. (|l|). This texture is purely planar, satisfies Hp — 0, but exhibits 
a vortex line running along the z-axis. The exchange energy density varies as 

Vex = ^ 10 

and diverges at r = 0. Therefore, we impose a short length cut-off of ro. Inside the vortex (r < rg) we assume a 
uniform energy density Ucore- The vortex radius tq can be chosen by balancing aM^/rg against Ucore at rg yielding 
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To ~ {aM"^ /UcoreY^'^ ■ Multiplying the core energy density by the cross-sectional area yields a vortex energy cost per 
unit length, C/y = Trr^ 

U core • 

This texture is probably not the absolute lowest in energy. We expect that the magnetization will rotate out-of- 
plane close to the vortex core. Such a distortion comes with a small cost in demagnetizing energy, but brings a large 
savings in exchange energy. Since this distortion is limited to the region close to the vortex p3, we may consider it 
as part of the core structure, and it is not relevant for the bulk droplet shape and energy. Large scale out-of-plane 
rotations cause large demagnetizing energy costs unless they are spatially nonuniform. Spatial nonuniformity will 
raise exchange energy costs. Hence, we expect the texture (|l|) is optimal sufficiently far from the vortex. 

The exchange energy can be reduced, and the vortex core energy eliminated altogether, by a change in topology 
from spherical to toroidal. Consider a torus with cylindrical symmetry about the z axis, with cross-section ±z(r) for 
Ri < r < Ro- Assuming texture (|^), and noting E^ore = E^^mag = 0, the total energy is the sum of surface and 
exchange energy 

/■^'^ , r^"" r^" r^'-''^ dz 

Etot = 2a del) rdr^JX + {dz/drf + — — / dcj) rdr —. (11) 

JO J Ri ^ Jo J Ri J -z(r) 

The curvature energy competes with the surface energy to determine the shape of the droplet in this model. For 
shapes with spherical topology, we must add an additional term to eq. (|ll| ) of the form UyLvi where Ly is the length 
of the vortex. However, since the exchange energy mimics the effect of the vortex core energy, we set [/y = in the 
following. 

The relative strength of each term is governed by a single dimensionless magnetization parameter defined so that 

= — — (12) 

(7 L 

where L is a measure of the linear droplet dimension (L {volumey-^^). Scaling lengths by L and energy by AttL'^u^ 
we introduce a dimensionless form of the total energy 

StoMr)]^ rdr^l + {dz/drf -f '-^ / rdr-, (13) 

J-Ri ^ JTli T 

where script quantities are dimensionless. 

To calculate the shape of the droplet we used the program Surface Evolver by Kenneth Brakke pO[ | . The program 
approximates a surface using a triangular mesh, and the vertices of the mesh evolve to minimize the total energy of 
the surface subject to various constraints of the problem (fixed volume in our case). 

We find that for small values of A4 the droplet takes a dimpled and slightly oblate "apple" shape (see Fig. |]). This 
occurs because the droplet reduces its curvature energy by reducing the length of the vortex (hence the dimpling 
where the vortex meets the surface) and by placing the bulk of the fluid as far as possible from vortex (hence the 
oblate shape). As M increases, the dimpling grows and the vortex shortens. This occurs even in the absence of a 
vortex core energy and is driven by the diverging exchange energy density near the vortex. For sufficiently large A4, 
the vortex shrinks to zero length and the droplet changes topology, acquiring a "donut" shape (see Fig. ||). 

For large M the hole at the center of the donut is large. As M decreases, the hole shrinks, and the walls of the hole 
become nearly vertical forming a cylinder of radius Ri. We can estimate this radius by balancing the Ri dependence 
of the cylinder surface energy against the Ri dependence of the magnetic energy outside the cylinder (and thus inside 
the torus). We estimate 

^.2,i. (14, 

and 

-dR- " -''^^^^ W ^ ^ 

from which it follows that Ri « aM^/2f7 (in dimensionless variables, TZi « ^A^^). Thus the donut hole remains for 
all values of A4, vanishing as ^ 0. Intriguingly, the structure of the donut hole for small A4 resembles the dimple 
on the apple surface at the same value of The dimple is a cylindrical hole of radius Ri, but the hole extends only 
part way through the apple core, with the vortex occupying the remainder. 

The energy of the apple is lower than the energy of the donut for small M and remains lower until just before the 
apple core shrinks to zero. For larger values of Ai, the apple shape does not exist and the donut is the stable shape. 
For a droplet with volume such that C = 12 and vortex radius ro = 0.1, the pinch-off point is close to M = 0.2. 
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IV. CONCLUSION 



In this paper, we fix a simple magnetization texture and then calculate the shape. Surface evolver is suited to 
optimization of shape. To express the total energy of the droplet as a function of its shape, we required that the 
magnetization be cylindrically symmetric and confined to a plane. The vortex core was assumed to be straight, and 
we neglected the energy of the vortex core. 

For a more rigorous study, a simultaneous calculation of the magnetization texture and the shape is needed. 
Since the magnetization may rotate out of plane near the vortex p^ , and possibly break the cylindrical symmetry, 
the analysis will require breaking up the volume of the droplet into finite elements and evolving the shape to 
minimize the sum of demagnetizing, curvature, vortex and surface energies. However, our simple analysis illustrates 
the nontrivial nature of the problem and gives an idea of shapes that might occur for ferromagnetic liquid droplets. 
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FIG. 1. Apple shape obtained by minimizing eq. ^ for = 0.17. 
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FIG. 2. Donut shape obtained by minimizing 



